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The  Effect  of  Cone  Angle  on  Penetration  Resistance 


Edward  A.  Nowatzki  and  Leslie  L.  Karafiath 


Research  Department 
Grumman  Aerospace  Corporation 
Bethpage,  New  York  11714 


ABSTRACT 

Static  and  dynamic  cone  penetration  tests  are  widely  used  in 
the  field  and  laboratory  to  obtain  comparative  measures  of  soil 
strength.  Theoretically  rigorous  methods  of  analysis  of  the  pene¬ 
tration  process  are  nonexistent;  approximate  methods  have  only 
limited  applicability. 

A  theoretically  correct  three  dimensional  analysis  of  cone 
penetration  using  plasticity  theory  and  the  Coulomb  failure  cri¬ 
terion  is  presented.  The  differential  equations  of  plastic  equi¬ 
librium  are  solved  numerically  for  an  ideal  uniform  dry  sand  to 
show  the  variation  of  slip  line  field  geometry  with  changes  in  the 
apex  angle  of  the  cone.  The  results  indicate  that  with  increasing 
apex  angle,  less  soil  volume  is  affected.  The  results  of  a  series 
of  laboratory  tests  using  variously  angled  cones  on  Jones  Beach 
sand  are  plotted  to  show  how  the  value  of  cone  index  increases 
with  increasing  apex  angle.  Cone  index  is  defined  as  the  ratio 
of  the  soil  resistance  to  the  base  area  of  the  cone.  The  differed 
tial  equations  of  plastic  equilibrium  are  again  solved  numerically 
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for  soil  and  boundary  conditions  that  correspond  to  those  of  the 
experiments.  These  results  support  the  validity  of  the  theory  in 
dense  sands.  The  experimental  results  also  demonstrate  that  soil 
compressibility  affects  the  oone  index  to  the  extent  that  it  no 
longer  serves  as  a  measure  of  frictional  strength.  For  loose 
soils,  differences  in  cone  angles  have  little  rffect  on  cone  in¬ 
dex,  all  other  conditions  being  equal.  To  identify  this  condition 
it  is  recommended  that  two  cones  be  used  in  penetrometer  tests, 
one  having  an  obtuse  apex  angle,  the  other  an  acute  apex  angle. 

The  theoretical  and  experimental  results  are  correlated  to 
show  how  the  theory  may  be  used  for  any  soil  to  predict  the  angle 
of  internal  friction. 

1 .  Introduction 

The  merits  of  using  penetrometer  data  for  determining  soil 
properties  have  been  discussed  extensively  [Fletcher,  G.  (1965)  . 
Meyerhof,  G.  G.  (1956),  Schultze,  E.  and  Knausenberger ,  H. 

(1957)].  In  many  of  these  reports  empirical  expressions  are  derived 
purporting  to  relate  penetration  test  parameters  to  soil  properties, 
for  example,  blow  count  data  of  the  Standard  Penetration  Test  (SPT) 
to  the  relative  density  of  the  soil  being  penetrated.  Although 
most  of  the  effort  has  been  directed  toward  dynamic  tests  such  as 
the  SPT,  some  consideration  has  also  been  given  to  relating  soil 
parameters  to  static  cone  penetration  characteristics  (WES,  1964). 
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At  any  rate,  very  little  attention  has  been  paid  to  a  theoretical 
analysis  of  the  interaction  between  a  cone  penetrometer  and  the 
failing  soil  during  the  penetration  process,  although  there  exists 
in  the  literature  a  well  established  basis  for  such  a  study.  The 
following  paragraph  presents  a  brief  review  of  the  pertinent  con¬ 
tributions  in  the  area  of  plasticity  analysis  of  soils*  problems . 

The  theory  of  static  equilibrium  has  been  combined  with  the 
Coulomb  failure  criterion  and  applied  to  studies  of  soli  bearing 
capacity.  Prandtl  (1920)  solved  the  resulting  differential  equa¬ 
tion  of  plastic  equilibrium  for  a  strip  footing  on  weightless  soil 
(dane  strain  condition) .  Cox,  Eason,  and  Hopkins  (1961)  developed 
a  general  theory  of  axially-symmetric  plastic  deformations  in  ideal 
soils  and  applied  it  to  the  problem  of  the  penetration  of  a  smooth, 
rigid,  flat-ended  circular  cylinder  into  a  semi-infinite  mass  of 
weightless  soil.  Drucker  and  Prager  (1952)  and  more  recently 
Spencer  (1962)  extended  the  theory  for  the  plane  strain  case  to 
include  body  forces  such  as  soil  weight  and  cohesion;  while  Cox 
(1962)  ar  .  more  recently  Larkin  (1968)  did  the  same  for  the  axially- 
symmetric  case.  In  cases  where  the  characteristic  relations  for 
the  governing  differential  equations  cannot  be  integrated  explicitly, 
numerical  methods  are  used.  Sokolovskii  (1965)  presented  the  most 
widely  used  procedure  of  numerical  integration,  a  finite  difference 
approximation  based  on  the  method  of  characteristics. 
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Studies  relating  specifically  to  cone  indentation  problems 
are  less  numerous.  Sneddon  (1965)  derived  a  solution  to  the  cone 
indentation  problem  within  the  framework  of  the  classical  theory 
of  elasticity.  Meyerhof  (1961)  and  Berezantzev  (1967)  offered 
approximate  solutions  of  the  axially-symmetric  problem.  The 
former  applied  the  solution  to  a  study  of  the  effect  of  surface 
roughness  on  cone  penetration;  the  latter  investigated  the  use  of 
penetrometer  data  to  determine  friction  angle. 

2 .  Objective 

An  analysis  is  presented  of  the  penetration  of  a  perfectly 
rigid  cone  into  an  ideal  granular  soil  whose  strength  properties 
are  defined  by  the  Coulomb  failure  criterion.  PTastic  stress  states 
are  considered  to  be  symmetric  with  the  central  axis  of  a  right 
circular  cone.  A  purely  frictional  soil  is  assumed  so  that 
local  pore  pressure  buildup  may  be  neglected.  This  assures  that 
failure  occurs  in  shear  zones  rather  than  along  a  single  failure 
surface.  In  the  model  used  in  the  analysis  the  limit  load,  ob¬ 
tained  from  a  solution  which  satisfies  the  basic  differential 
equations  of  plasticity  and  the  boundary  stress  condition,  is  con¬ 
sidered  a  lower  bound  (Drucker  and  Prager,  1952) . 
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3 .  Scope 

This  study  is  arranged  in  the  following  order.  In  Section  4 
the  variables  are  identified  and  defined.  Section  5  contains  the 
governing  equations  of  plastic  equilibrium  for  the  axially-symrueLiic 
case.  In  Section  6  the  limitations  of  the  present  investigation 
are  discussed.  The  major  results  obtained  from  the  theory  within 
those  limitations  are  presented  and  discussed  with  reference  to 
experimental  data.  Conclusions  are  drawn  at  this  point  in  the  analy¬ 
sis.  Finally,  in  Section  7  the  results  of  this  study  are  summarized 
and  remarks  made  concerning  them  and  their  relation  to  future  research 
in  this  area. 

4 .  Notation 

The  quantities  defining  the  geometry  of  the  problem  are  given 
in  Fig.  1. 

A  -  area  of  cone  base 
o 

c  ■  cohesion 

G  *  slip  line  geometry  similitude  factor 

Cl  «  cone  index,  defined  in  psi  as  the  penetration 
resistance/AQ  where,  in  this  study,  the  resis¬ 
tance  at  6-in.  penetration  is  used  as  a 
reference 

R  ■  radius  of  cone  base 
o 

w  ■  surcharge 
r,z  ■  coordinates 
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»  depth  to  which  base  of  c  ne  has  penetrated 
a  ■  apex  angle  of  cone  or  cot  e  angle 
0  *  complement  to  apex  semi-fngln 

y  *=  unit  weight  of  soil 
6  ■  friction  angle  between  cone  a?.d  soil 

6  «  angle  between  r  axis  and  r:ajor  principal  stress 

cp  *=  angle  of  internal  friction  ot  soil 

Ip  ■*  c  cot  cp 

p.  =  tt/4  -  cp/2 

a  «  (a^  +  a.j)/2  +  ip  (in  general) 

5 .  Formulation 
* 

The  following  set  of  differential  equations  represents  Che 
theoretically  rigorous  formulation  to  the  probltn.  of  determining 
axially-symmetric  plastic  stress  states  and  slip  r,e  fields. 

da  +  2 a  tan  <t>d©  -  (^sin(+  <t>)dr  +  cos(+  cp)dzy  ■» 

--  (^sin  4>dr  +  tan  <t>(l  -  sin  4>)dzj  «  0 


dz  -  dr  tan(0  +  p)  ^1) 

This  set  of  equations  is  obtained  by  combining  the  equations  of 
equilibrium  derived  from  plasticity  theory  with  the  Coulomb  fail¬ 
ure  criterion.  The  circumferential  stress  is  assumed  to  be  the 
intermediate  principal  stress  and  to  be  equal  to  the  minor  principal 
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stress  (o^).  The  r-axis  of  the  coordinate  system  is  parallel 
to  the  ground  surface  and  the  positive  z_axis  is  perpendicular 
to  it  into  the  soil  mass  (see  Fig.  1).  The  upper  sign  refers  to 
the  family  of  slip  lines  corresponding  to  the  first  characteristics 
of  the  differential  equations  (i-lines),  and  the  lower  sign  to 
the  second  (j-lines). 

No  closed  form  solution  to  these  equations  exists.  Several 
numerical  solutions  have  been  presented;  however,  these  are  re¬ 
stricted  to  the  axialiy-symmetric  surface  loading  of  the  semi¬ 
infinite  half  space.  The  equations  below  are  the  numerical  form 
of  Eqs ,  (1)  used  to  study  the  effect  of  soil  properties  and  cone 
parameters  on  the  penetration  characteristics  of  a  right  circular 
cone  p  netrating  soil.  To  keep  the  problem  perfectly  general, 
soil  body  forces  have  been  included  in  the  formulation.  For  a 
given  set  of  loading  conditions  (w)  over  the  horizontal  soil 

surface,  the  values  of  r.  .,  z.  .,  a.  .,  and  9.  .  are  computed 

i » J  1 ,  J  i » J  1 » J 

for  an  adjacent  nodal  point  (slip  line  intersection  point;  refer  to 
Fig.  1)  by  use  of  the  set  of  recurrence  relationships: 


i  j  "*  (zi_i  j  "  z ^  j“l  6 lri  j-1  ”  e2ri-l  j^^^^l  ™  (  2 ) 


z .  .  -  z .  .  .  +  (r .  .  -  r,  ,  .  )e0 
i,J  i-l, J  i,J  i-l.J  2 


(3) 


where  r.  .  and  z.  .  are  the  coordinates  of  the  adjacent 

i,j  i,j 

nodal  point,  and  =  tan(ti^  +  p),  €2  *  tan(0^_^  ^  -  p) .  With 
these  values  of  r.  .  and  z.  .  the  computation  is  continued  for 
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ai,j  "  +  tan  ’  6i-l,j))  + 

+  °i,j-lC  "  +  A/ri-l,j)]  / 


(ai,j-l  + 


,  .  ■  a.  .  ,  -  a.  ,  ,  +  2  tan  qp(c.  .  .  +  a.  ,  ,  . 

i,j  i,j'l  i'l.j  i»J-l  i.j'l  i‘l>J  i"l,J 

+  D  -  C  +  '  ai,j-lB/'ri,j.lj  / 


2  tan  +  VlJ* 


where 


A  -  sin  cp  (r  ±  ^  -  r^  j)  “  tan  ‘X1  *  sin  <P)  (Zj.  j  -  z^  j) 

B  -  sin  9(ri^  -r^^)  +  Can  ^ C 1  -  sin  cpXz^  -  z^^) 

c  -  Hzi,j  -  - tan  -  ri-i,J)) 

D  -  Hzi,j  -  zi,j.i +  can  '"('i.j  ■  ri,j-i>) 


To  apply  these  recurrence  relationships  to  the  problem  of 
cone  indentation,  the  geometric  boundary  conditions  as  well  as 
the  stress  boundary  conditions  had  to  be  formulated  appropriately 
and  included  in  the  computer  program.  The  geometric  boundary  con¬ 
ditions  are  simply  mathematical  descriptions  of  the  cone  geometry 
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and  its  position  at  depth  in  terms  of  r,  z,  a,  and  Rq.  The 
stress  boundary  conditions  and  the  method  of  computation  are 
essentially  the  same  as  those  described  by  the  authors  elsewhere 
for  two  dimensional  conditions  (Karafiath  and  Nowatzki,  1970). 

Briefly,  the  stress  boundary  condition  on  the  horizontal  plane 
through  the  base  of  the  cone  is  given  by  the  surcharge  (w)  and 
the  overburden  soil  pressure  (see  Fig.  1).  Overburden  shear  is 
disregarded.  The  slip-line  field  in  the  passive  zones  is  computed 
by  Eqs .  (1)  starting  with  these  boundary  values  and  an  assumed 
value  for  the  horizontal  extent  of  the  passive  zone.  In  the 
radial  shear  zone,  the  same  equations  are  used,  but  special  con¬ 
sideration  is  given  to  the  central  point  where  the  j -lines  con¬ 
verge  (point  Q  in  Fig.  1).  This  point  is  a  degenerated  slip 
line,  where  8  changes  from  the  value  at  the  passive  boundary  to 
that  specified  at  the  active  zone  boundary.  The  total  change  in 
0  is  divided  by  the  number  of  slip  lines  converging  at  this  point 
to  obtain  an  equal  A 0  increment  between  two  adjacent  slip  lines. 
The  o  values  for  each  increment  are  computed  from  the  equation 

o  =  o  Po)tan  ^ ,  which  is  the  solution  to  Eqs .  (1)  if  both 

o 

dr  and  dz  vanish.  Allowance  must  be  made  for  the  angle  3  as 
well  as  6  in  assigning  transition  values  of  8  and  o  between 
the  active  and  passive  zones.  With  these  values  of  8  and  a  for 
each  slip  line  at  this  point,  the  coordinates  as  well  as  the  a 
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and  6  values  for  all  other  points  in  the  radial  shear  zone  can 

be  computed  by  Eqs.  (1).  In  the  active  zone  the  same  equations 

are  used  except  for  the  points  at  the  loaded  surface  of  the  cone 

itself,  where  8.  ,  is  assigned  and  z.  .  -  z  +  (R  -  r,  .)  tan  u. 

1  >  J  °  °  1  >  J 

Here: 


i  >  j 


•  t  .  tan(  (0,  ,+  8  .  . )  /  2  -  p. }  -  z .  ,  .  +  R  tan  a  +  z 

i-1  llAiA _ _ J. _ liLj _ 2 _ 2 


tan 


(<ei,j  +  9i-ij)/2 


p.  j  +  tan  a 


(10) 


'l,j  -  °i-lj  +  tan  9  +  c  •  r 


o,  ,  .A 


(11) 


The  numerical  computation  is  performed  and  adjustments  made,  if 
necessary,  to  the  value  assumed  for  the  horizontal  extent  of  the 
passive  zone  until  the  slip  line  field  "closes"  on  the  axis  of 
symmetry  at  the  apex  of  the  cone. 

6  .  Results  and  Conclusions 

a.  Theoretical  Results  -  Slip  Line  Field  Geometries 
The  numerical  computation  of  the  slip  line  field  geometries 
and  associated  stresses  by  the  recurrence  relationships  [Eqs.  (2) 
through  (5)]  was  performed  on  an  Adage  Inc.  time “sharing  computer 
system.  This  system  is  based  on  Digital  Equipment  Corporation 
PDP-10  processors.  To  show  specifically  the  effect  of  cone  angle 
on  the  slip  line  field  geometry,  the  governing  differential  equations 
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were  solved  for  a  set  of  ideal  soil  conditions  that  describe  a 
homogeneous,  dry,  and  purely  frictional  sand.  These  conditions 
are  c  «  0,  <p  -  37°,  y  -  100  pcf,  w  -  1  psf,  and  6  ■  20°.  The 
numerical  results  were  plotted  automatically  and  electronically 
on  the  display  tube  of  a  Computer  Displays  Inc.  Advanced  Remote 
Display  System.  The  slip  line  fields  for  Rq  -  0.034  ft  (radius 
of  the  WES  cone)  and  a  -  15.5°,  30°,  60°,  90°,  120°,  and  150° 
have  been  reproduced  and  are  shown  in  Figs.  2.  Since  the  problem 
is  axially-symmetric, only  half  the  total  field  has  been  shown; 
the  dashed  line  indicates  the  central  axis  of  the  cone.  The 
geometric  scale  may  be  obtained  for  each  figure  from  the  knowledge 
that  the  base  radius  of  the  cone  is  0.4  in.  The  scale  in  Fig.  2a 
is  four  times  that  in  Fig.  2b. 

For  given  soil  strength  parameters  c  and  cp  the  slip  line 
fields  representing  the  solution  of  the  differential  equations  are 
geometrically  similar  only  if  the  ratio  G  =  yll  (c+w  tan  cp)  is 
the  same  (Cox,  1962),  where  l  is  a  characteristic  length  usually 
taken  equal  to  Rq.  Although  the  slip  line  geometry  similitude 
factor  G.  as  defined  by  Cox  for  the  case  of  bearing  capacity,  is 
the  same  (4.512)  in  all  of  the  cases  represented  in  Figs.  2,  it 
is  obvious  that  the  slip  line  geometries  differ,  an  .  ndication 
that  the  G-equality  is  a  necessary  but  not  sufficient  condition 
for  slip  line  field  similitude.  In  addition  to  G,  there  must 
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also  be  equality  In  6  and  a  similitude  In  the  geometric  condi¬ 
tion  at  both  the  free  and  loaded  boundaries  to  obtain  complete 
geometric  similitude  of  the  slip  line  field. 

The  salient  feature  of  Figp.  2  is  that  they  distinctly  show  a 
contraction  of  the  radial  shear  zone  with  decrease  in  cone  angle. 

In  Fig.  2b  the  contraction  is  so  pronounced  that  the  individual 
i-  and  j-lines  are  hardly  discernible  with  the  scale  used.  The 
active,  passive,  and  radial  shear  zones  all  have  curvilinear 
boundaries  due  to  the  three  dimensional  nature  of  the  problem,  an 
indication  that  the  geometries  obtained  from  the  solution  of  the 
theoretically-correct  differential  equations  differ  from  those 
obtained  by  using  the  Prandtl  solution  foi;  weightless  soil  and  the 
log-spiral  approximation  in  the  radial  zone. 

Also  indicated  by  Figs.  2  is  the  fact  that,  with  decreasing 
cone  angle,  the  affected  volume  of  soil  increases.  For  the  soil 
and  cones  used  in  Figs.  2,  the  volume  of  the  body  of  revolution 
formed  by  the  slip  line  field  for  a  =  30°  is  about  ten  times 
greater  than  that  for  cr  =  150°.  For  compressible  soils,  the  size 
of  the  affected  mass  directly  influences  the  load-penetration  re¬ 
lationship.  The  material  must  be  compressed  to  a  state  in  which 
friction  is  fully  mobilized  before  shear  failure  along  slip  lines 
can  take  place.  The  load  necessary  to  accomplish  this  compression 
is  usually  less  than  the  load  needed  to  fail  the  soil  in  shear.  Ic 
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follows  that  the  larger  the  volume  of  the  slip  line  field,  the 
more  the  soil  mass  must  be  compressed  to  mobilize  the  friction 
fully.  Therefore,  cone  indices  obtained  with  cone  shapes  that 
result  in  large-volume  slip  line  fields  are  likely  to  be  less 
representative  of  the  Coulombic  strength  than  of  the  compressi¬ 
bility  of  the  material. 

b.  Experimental  Results 

A  series  of  penetration  tests  were  conducted  on  Jones  Beach 

sand  using  aluminum  cones.  The  friction  angle  between  the  aluminum 

surface  of  the  cone  and  the  sand  was  assumed  to  be  uniform  and  equal 

to  15  degrees  on  the  basis  of  experiments  performed  by  Mohr  and 

Karafiath  (1967)  .  The  purpose  of  the  tests  was  to  determine  the 

effect  of  cone  angle  (a)  on  the  value  of  the  cone  index  (Cl)  . 

The  base  area  of  cones  having  apex  angles  of  150,  90,  and  30  de- 

2 

grees  was  0.-5~in.  .  The  30“degree  cone  corresponds  to  the  WES 

2 

cone.  Also  used  was  a  cone  having  a  base  area  of  1.04-in.  and 
an  apex  angle  of  15.5  degrees.  This  cone  corresponds  to  that  used 
until  1956  by  the  North  Dakota  State  Highway  Department  for  flexi¬ 
ble  pavement  design.  The  physical  properties  of  Jones  Beach  sand 
are  listed  in  Table  1.  Uniform  soil  beds  were  prepared  to  a  nar¬ 
row  range  of  desired  densities  in  a  facility  specially  designed  for 


this  purpose . 
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To  conduct  the  tests,  the  cones  were  attached  interchangeably 
to  a  12-in.  rod  and  the  entire  assembly  mounted  on  the  loading 
frame  of  an  Instron  testing  instrument  Model  TM-M  (Fig.  3).  The 
rate  of  penetration  was  set  at  10  cm/min.  Load-penetration  curves 
were  obtained  automatically  on  a  synchronized  strip-chart  recorder. 
Values  of  Cl  as  defined  for  this  study  were  determined  directly 
from  the  load-penetration  curves.  The  results  of  the  penetration 
test  series  are  given  in  Figs.  4  and 

Figure  4  shows  the  variation  of  cone  index  with  change  in 
cone  angle  for  Jones  Beach  sand  at  various  relative  densities. 

It  is  clear  from  the  figure  that  v>hen  the  material  is  loose 
(88  pcf  <  /  <  94  pcf)  the  cone  angle  has  very  little  effect  on 
the  cone  index  (all  values  of  Cl  are  less  than  approximately  9 
regardless  of  the  size  of  the  cone  angle).  On  the  other  hand,  for 
dense  materials  (100  <  y  <  105)  Cl  varies  from  approximately  19 
for  the  15.5  degree  cone  to  between  44  and  60  for  the  150  de¬ 
gree  cone.  These  results  seem  to  verify  the  previously  discussed 
effect  of  compressibility  on  the  value  of  cone  index.  Apparently 
the  frictional  strength  of  the  loose  material  cannot  be  fully 
mobilized  until  the  material  is  sufficiently  compressed  to  allow 
for  complete  shear  failure.  Calculations  made  on  the  basis  of 
volune  change  properties  of  the  Jones  Beach  sand  indicated  that 
the  average  percent  volume  change  of  the  soil  mass  within  the  slip 
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line  field  was  virtually  independent  of  the  apex  angle  of  cones 
having  the  same  base  area.  Consequently,  the  displacement  neces¬ 
sary  to  mobilize  the  full  friction  is  roughly  proportional  to  the 
volume  of  the  slip  line  field.  These  results  suggest  that, in 
order  to  use  the  cone  index  as  a  valid  measure  of  the  frictional 
strength  of  a  soil,  two  cones  should  be  used,  both  having  the 
same  base  area  but  one  having  an  obtuse  apex  angle,  the  other  an 
acute  apex  angle.  If  the  values  of  Cl  obtained  from  these  two 
cones  are  similar,  then  the  penetration  resistance  is  governed  by 
the  compressibility  of  the  soil.  If  on  the  other  hand,  the  two 
values  are  markedly  different,  then  a  relationship  may  properly  be 
sought  between  Cl  and  the  strength  parameters  of  the  material. 

Of  course,  the  problem  is  not  entirely  that  simple.  In  addi¬ 
tion  to  the  coup  ressibility  effect  at  lower  relative  densities, 
there  is  the  effect  of  the  variation  of  <p  with  y.  Figure  5 
shows  for  each  of  the  cones  investigated  the  change  in  cone  index 
with  variation  of  unit  weight.  It  is  impossible  from  such  a  plot 
to  distinguish  which  of  the  two  effects  has  the  greater  influence. 
However,  it  seems  reasonable  to  assume  that,  in  the  range  of  unit 
weights  over  which  the  material  may  be  considered  relatively  dense, 
the  compressional  effect  is  negligible. 

c.  Comparison  between  Experimental  and  Theoretical  Results 

In  Fig.  6  the  experimentally  obtained  values  of  Cl  are  com¬ 
pared  over  a  wide  range  of  unit  weights  with  the  values  determined 
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theoretically  by  the  solution  of  the  differential  equations  of 
plastic  equilibrium.  For  clarity,  only  the  results  for  the  150  de¬ 
gree  cone  are  presented.  The  theoretical  curves  suggest  that  for 
<p  constant  there  is  little  change  in  Cl  with  variation  in  7. 

The  theoretical  curves  also  show  that  for  a  given  unit  weight  the 
value  of  Cl  is  very  sensitive  to  changes  in  <p.  The  experimental 
curve,  which  shows  a  pronounced  decrease  of  Cl  with  decrease  in 
7,  not  only  reflects  the  change  in  <p  with  unit  weight  and  stress 
level  (refer  to  Table  1),  but  also  includes  the  effect  of  soil  com¬ 
pressibility  discussed  above.  Unfortunately,  the  experimental 
curve  itself  does  not  distinguish  between  these  effects;  however, 
for  reasons  cited  above  it  seems  that  the  compressibility  has  the 
least  influence  on  tests  performed  with  dense  material.  For  ex¬ 
ample,  from  Fig.  6  a  cone  index  of  43  for  Jones  Beach  sand  at 
103  pcf  indicates  a  cp-angle  of  approximately  38-39  degrees. 

For  the  stress  level  involved,  this  value  of  cp  determined  from 
cone  penetrometer  data  agrees  quite  well  with  the  values  obtained 
from  triaxial  tests  listed  in  Table  1. 

Therefore,  although  the  theory  of  plastic  equilibrium  has  not 
been  modified  in  this  study  to  include  in  the  computation  of  Cl 
the  effects  of  a  curvilinear  Mohr  envelope  (Berezantzev  and  Kovalev, 
1968;  Szymanski,  1958),  or  the  contribution  of  soil  compressibility, 
we  believe  that  curves  such  as  those  presented  in  Fig.  6  can  be  used 
to  estimate  an  average  value  of  q?  from  cone  penetrometer  data. 
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7 .  Summary  and  Discussion 

The  major  results  of  the  present  Investigation  are  summarized 
as  follows: 

a.  For  the  axia lly- symmetric  case,  the  slip  line  field  geome¬ 
try  derived  from  the  theory  of  plastic  equilibrium  for  a  dry,  uni¬ 
form  sand  being  penetrated  by  a  cone  differs  markedly  from  that 
obtained  by  the  Prandtl  solution  for  weightless  soil  and  the  log- 
spiral  approximation. 

b.  For  materials  at  high  relative  densities,  the  cone  index 
varies  significantly  with  the  size  of  the  penetrometer  apex  angle, 
all  other  conditions  being  equal.  This  is  not  observed  for  loose 
materials.  Soil  compressibility,  the  variation  of  cp  with  y, 
and  the  curvilinear  nature  of  the  Mohr  failure  envelope  can  account 
for  this  difference  in  performance. 

c.  Solutions  derived  from  the  theory  of  plastic  equilibrium 
agree  well  with  experimentally  obtained  cone  penetrometer  data  and 
can  be  used,  under  certain  conditions,  to  estimate  the  strength 
parameter  cp  for  a  dense  dry  sand.  In  all  penetrometer  investiga¬ 
tions  the  use  of  two  cones  is  recommended  to  avoid  misinterpreta¬ 
tion  of  the  cone  index.  One  cone  should  have  an  obtuse  apex  angle, 
the  other  an  acute  angle. 
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Further  investigations  are  needed  to  enhance  the  theory  pre¬ 
sented  in  this  study.  For  example,  it  would  be  very  desirable  to 
incorporate  into  the  computer  program  the  curvi linearity  of  the 
Mohr  envelope  and  the  dependence  of  <p  on  unit  weight.  Similarly, 
the  effect  of  compressibility  on  the  penetration  resistance  of  a 
material  should  be  expressed  quantitatively,  and  criteria  for  the 
mobilization  of  friction  in  the  slip  line  field  established. 
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TABLE  1 

ANGLE  OF  INTERNAL  FRICTION  OBTAINED  FROM 
TRIAXIAL  TESTS  ON  JONES  BEACH  SAND 


Range  of  9 

Range  of  y  over  Normal  Stress  Levels  of  Interest 


103.5 

-  105.5  pcf 

41°  -  37° 

95.5 

-  97.5  pcf 

38°  -  30° 
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?ig-  Geometry  for  Cone  Indentation  Problem.  Boundary  conditions  an 
slip  lines  are  symmetric  with  the  z-axis. 
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Fig.  2  Slip  Line  Fields  for  Cones  Penetrating  Ideal  Soil 
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Fig.  4  Cone  Index  at  6 -in.  Penetration  versus  Cone  Angle 

for  Jones  Beach  Sand  at  Various  Relative  Densities. 
Experimental  results. 


Fig.  6  Cone  Index  at  6-Ln.  Penetration  versus  Unit  Weight  for 
150°  Cone.  Experimental  and  theoretical  results. 


